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Abstract 

We consider the effective Hamiltonian of four quark operators in the Standard 
Model in the exclusive and quasi-inclusive decays of the type B — ► K^rf, B — ► n'X s , 
where X s contains a single Kaon. Working in the factorization assumption we find 
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that the four quark operators can account for the recently measured exclusive decays 
B — ► rj'(rj)K and B — ► i^7r for appropriate choice of form factors but cannot explain 
the large quasi-inclusive rate. 

1 Introduction 

Recently CLEO has reported a branching ratio for the process B~ -> t/-K" _ at (7.ll|f ± 
0.9) x 1CT 5 IQ. Upper limits on related exclusive decay modes with 77 or 77 with a X* 
mesons in the final states were also obtained. The quasi-inclusive process B — > i]'X s for high 
momentum 77' was also measured with a high rate @ 

B^ri'X s = (6.2 ± 1.6 ± 1.3) x l(T 4 0y > 2.0GeV) (1) 

where X s stands for one Kaon and up to 4 pions with at most one 7r°. 

Several explanations for the large quasi-inclusive decay rate have been proposed both 
within the Standard Model and beyond []3], f|, |5|, However, some of these analyses have 
underestimated the effect of the effective four quark operators to this process. For example, 
it has been assumed that this contribution is small based on the smallness of Kj&B- However, 
the penguin contribution to this process cannot be neglected. A simple estimate based on 
the magnitudes of the CKM elements associated with the tree and the penguin diagrams 
clearly suggests that this process is dominated by penguin contributions if one considers 
only the four quark effective Hamiltonian. In this paper we calculate the contribution of the 
effective Hamiltonian of the four quark operators to the exclusive decays B — > K^*'rj'{r]) and 
B — > Kir along with the quasi-inclusive decay of the B — > rj'X s . 

In the sections which follow, we describe the effective Hamiltonian of four quark operators, 
our calculation of the the exclusive and quasi-inclusive rates. 
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2 Effective Hamiltonian 



In the Standard Model (SM) the amplitudes for hadronic B decays of the type b — > qff are 
generated by the following effective Hamiltonian @: 
/-» 10 

H h = -^WfbVf^cxOif + c 2 oi f ) - + v cb v; q ct + v tb v*4)o«) + h.c. (2) 

V Z i=3 

where the superscript u, c, t indicates the internal quark, / can be u or c quark and q can 
be either a d or a s quark depending on whether the decay is a AS = or AS = — 1 process. 
The operators 0\ are defined as 

0% = WpLffiffrfLb* O q 2f = q llM LffrLb 

°h = q r r» L bq' r rv L ( R )<l' °l,e = Qal^Lb^f.LiR)^ (3) 

3 3 
Or, 9 = -q-f^Lbe q/ qYR(L)q Of >10 = -q a ^Lbpe q/ q^R(L)q a 

where R(L) = 1 ± 75, and q' is summed over u, d, and s. 0\ are the tree level and QCD 
corrected operators. are the strong gluon induced penguin operators, and operators 

O 7 _i are due to 7 and Z exchange (electroweak penguins), and "box" diagrams at loop 
level. The Wilson coefficients c{ are defined at the scale \i ~ and have been evaluated 
to next-to- leading order in QCD. The c\ are the regularization scheme independent values 
obtained in Ref. ||. We give the non-zero c{ below for m t = 176 GeV, a s {m z ) = 0.117, and 
/1 = m h = 5 GeV, 

d = -0.307 , c 2 = 1.147 , 4 = °- 017 > 4 = -0.037 , c\ = 0.010 , c\ = -0.045 , 
4 = -1.24 x 10" 5 , 4 = 3.77 x 10~ 4 , c\ = -0.010 , c\ = 2.06 x 10~ 3 , 

4; 5 C = -<$$/Nc = pr/*o 4; 9 c = p;' c , 4 ;i c = o (4) 

where N c is the number of color. The leading contributions to P l se are given by: P l s = 
(§)c 2 (f+GK,// )g 2 ))andP e i = (^)(jV cCl +c 2 )(f +G(m i)M 2 )). The function G(m, n, q 2 ) 
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is given by 



G(m, fi, q ) = A / x(l-x)ln — — dx (5) 

Jo /i 2 



All the above coefficients are obtained up to one loop order in electroweak interactions. The 
momentum q is the momentum carried by the virtual gluon in the penguin diagram. When 
q 2 > 4m 2 , G(m, fi,q 2 ) develops an imaginary part. In our calculation, we use m u = 5 MeV, 
m d = 7 MeV, m s = 200 MeV, m c = 1.35 GeV g [H| and use q 2 = m 2 /2. 



3 Matrix Elements for B K^f] f (rj), B —> Kit and B 

The effective Hamiltonian described in the previous section consists of operators with a 
current x current structure. Pairs of such operators can be expressed in terms of color singlet 
and color octet structures which lead to color singlet and color octet matrix elements. We 
use the factorization approximation, where one separates out the currents in the operators by 
inserting the vacuum state and neglecting any QCD interactions between the two currents. 
The basis for this approximation is that, if the quark pair created by one of the currents 
carries large energy then it will not have significant QCD interactions. Factorization appears 



to describe nonleptonic B decays rather well|TT|] . To accommodate some deviation from this 



approximation one can treat N c , the number of colors that enter in the calculation of the 
matrix elements, as a free parameter though the value of N c ~ 2 is suggested by experimental 
data on low multiplicity hadronic B decays. In this section we describe the calculation of 
matrix elements for the exclusive decays B — > rj'(rj)K^*\ B — > Kir and the quasi-inclusive 
decay B — > r)'X s . 
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V 




cos 9 — sin 9 




Vs 


rf _ 




sin 9 cos 9 







The 7] and rj' mesons are mixtures of singlet and octet states rji and r/g of SU(3). 

(6) 

uu + dd — 2ss (7) 
uu + dd + ss (8) 

where the mixing angle 9 lies between —10° and — 20°fl(|. We can express the decay constants 



" 8 = 7? 



fn'y ft 1 -, fv> m terms of octet and singlet decay constants fi,fs 



eu rd 

J r}' J rf 

fr)' 



3 u 

3^ 



fi sin 9 f 8 
cos 9— + —prj- 

Jit V £ In 



cos#y- — \/2 sin 9^j- 

Jn Jn 



(9) 
(10) 



and similarly one has 



/u rd 

Jr) 



cos 9 



■ a^ 1 
sin 9— 



V2 ' fn 

V2 cos Q-£- + sin 9^- 



11) 



'12) 



In the SU(3) and Nonet symmetry limit, the relations fs — fn — 130 MeV and j\ = hold. 
However these symmetries are not exact and we will consider values for fi and/s away from 
this symmetry limit. For example the value of fx ~ 1-28/^ indicates about a 30% SU(3) 
breaking effect. 

3.1 Exclusive Decay 

In this section we calculate the rates for the two body noleptonic decay for B — > r]'(r])K and 
B -> tcK . 

To evaluate the rates for the exclusive decays we have to calculate matrix element of the 
type < hK\H e ff\B > where h is a 7r or 7/(77) and H e ff, as already described in the previous 
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section, has the form. 

H eff = V u { Cl Oi + c 2 2 ) - ]T [v u c? + V£ + V t c\ 

i=3 

with V u = V* s V ub , V c = V*V cb and V t = V*V tb . 



(13) 



As an input to our calculation we need the form factors defined through 



< K(p K )\s^(l - j 5 )b\B{p B ) > 



m n — m 



+ 



(pb + Pk)v - 



K 



-?mW) 



(14) 



where q = p B — Pk- 

The matrix element < h\u^^{\ — 7s)fe|-B > is similarly described in terms of form factors 
F/V) and Ffr(q 2 ). 

For the K* modes we define 



.A* 



(15) 



< ir A |s7 M (l - 75 )m|0 >= m K *g K *e* 
The form factors in this case are defined as 

<K*(p')\r\B(p)> = W f ] e^elp a p' p + i{m B + m^A^k 2 ^ - is* ■ k ' M "~~ 



ttl b + m K * 
+ fe . h 2m K ,(A (k 2 ) - A 3 (k 2 )) kU 
k 2 



m B + m K * 



(p+pT 

(16) 



with 



A 3 (k 2 ) = 



(m B + m K *)Ai - (m B - m K *)A 2 
2m K * 



where k = p B — Pk*- 



3.1.1 B ± -> K^'irj) 



The amplitude for B ^ K rj' can be written as 
V2 



M 



V u (a iri Q K + a 2 Q v ,)- ]T ^{( T i>i + T % 2 r 2 )Q K + T*Q V ,} 



u,c,t 



(17) 



where 



t\ = 24 -24 --4 + -4 



™^ 1 , 1 , 1 



T~ = al + al — al + (2al - al) ; — + -al al a l m 

2 3 4 5 v e SJ 2m s ( y m b - m s ) 2 7 2 9 2 10 

T* = 4 + 2(4 + 4)^ + a\ 

m u + m s m b - m u 

with ai = Cl + f , a 2 = c 2 + §- c , 4 = c) + 4 +1 = c) +l + jf c , n = f r 2 = f , 
Q K = iFf '{ml,){m 2 B - m 2 K )f n , Q rf = iF$\m 2 K ){m 2 B - m 2 ,)f K , V { = V u ,V c ,V t and N c is 
effective number of colors. 

In the above equations we have used the quark equations of motion to simplify certain 
matrix elements. The expression for the amplitude can also be used for B — > rjK by making 
the necessary changes. It is also straight forward to write down the amplitudes for B — > K*rj' 
and B — > K*rj decays. 



M 
where 



G p 

V2 



V u (a 1 ft,A + a 2 m K *g K *B) - £ V t {(7* ft + T 2 ft)A + T 3 m K .g K .B} 



i=u,c,t 



(18) 



1 . 1 . 
24 - 24 - -4 + -4 



n = 4 +4 -4 -(24 



m 2 , 



2m s (mt ) + m s ) 2 



1 i 1 i 1 i 

^~ r>^7 0^9 o^lO 



T* = 4 + a 
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with A = 2m K *A e* ■ ps, B = 2e* ■ VbF\ (m 2 K*) an d we wm use 9 k* = 221 MeV. A 
similar expression for B — > can also be obtained. 



3.1.2 B ^ Kir 



We give below expressions for the various B — > i^7r decays 
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M(B~ -> K-k' 



V2 K 



i o a 10 + 



2m\ 



-(al - -al) 

{m s + m d )(m b - m d ) 2 



Q, 



if K (m 2 B -mi)F^(mj,) 



M{B- -> /TV 

M(i? -> 7T + ) 

Qtt 



G p 

71 



V5 



1 , 
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2m K *f K *e*.p B F^{m 2 K 



V2 



V u a 2 - \ \ <i\ r/, 

2 



2m^ 



^ (m s + m u )(m b - m u ) 

G F 



a fi + do 



M(5° -> K-*7T+) = -^[y„a 2 -V r i {oj + ai }]2m if ./jc.e*.pBFi B " +,r (m^ 



(19) 

Qtt 

(20) 



-7T ) 



M(S~ -> iT~7r° 



Gp 

V2 
G f 

V2 



V u {a x Q K + a 2 Q n } - y {t*Q k + T^Q n } 
V u {aiQ K > + a 2 Q n ,} - V t {t(Q k , + if Q v } 



(21) 



where 7? = -fa* + fa* = 7* Tj = 4 + < + ( ms+ jJ^) (4 + 4). T 2 = < + a\ , Q n 



tf K (m 2 B -ml)F B ^°(m 2 K ) and Q x = i^(m\-m\)F^ K {m\) Q„, = 2m K * f K *s* .p B F{ 
and Q K > = 2m K *e*.p B A^ K *^ 

3.2 Quasi-inclusive Decay 



B- 



The technique to handle quasi-inclusive decays have been described in detail in Ref . [|12| . We 
represent the total amplitude as the sum of a "two body" and a "three body" piece. The 
"two body" amplitude is given by 

. Gp 



Mo 



:-|#p2 < X\*f{l - ls)b\B > FL U 
:|< M -7a)^>g<, 



(22) 



where 



FL, 



K|ci + |) + ^{ 2+ ( 1 + i_) r } + A (i. + (_L + llr 
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4{ i -^ +i )h + HiH^ +1 M- 4i ° 

with Ai = -{ti}V u + ciV c + cjV t ) and r = ^ 
The "three body" piece has the form 

M 3 = <r ] \uY(l--f 5 )b\B><X\s-f fl (l-^)u\0> L u 

+ < T)\u{l - 7s)fe|5 >< X\s(l + 7s)m|0 > £ u (23) 



with 



iV c V ViV c V V.V, 
Details of the calculations of the branching fraction and decay distributions are given in Ref 







4 Results and Discussion 

For the exclusive decays to the rj' the inputs to the calculation are the form factors, the values 
of f5 or alternately the values of f\ and f$, N c and the mixing angle. We will try to fit 
the experimental number for B —>■ Krj' by assuming fi ~ (1.0 — 1.5)/^ and /g ~ (1.0 — 1.5)/,,-. 
We will take N c = 2, the effective number of colors. We find that there are several solutions 
corresponding to different values for the set of parameters fx, f 8,0 that can reproduce the 
experimental data on B — > Kr\ and the upper limit of B — > Kt] if we use the form factors 

in Ref [|l3|l For the CKM parameters we choose two sets (p = 0.15 and rj = 0.33) and 

4 For the BSW model the form factor Fq(0) for B — > rj and B — > rj are approximately same and so we 
assume this equality for the form factors in Ref [|l3| where only the B — > r) form factor is calculated. 
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(p = —0.15 and 77 = 0.33) JT^|. In Table. 1 we give the rates for the decays involving a 
and 77, rf and n in the final state for 6 = —17° and N c = 2 for two sets of the CKM 
parameters with various values of f\ and / 8 . Phenomenological studies involving radiative 
decays of the rj and rf indicate values for fi ~ l.l/ w and / 8 ~ 1.3/„- JTB|. We show the results 
of using the form factors in Ref [13] and Ref 0] in the third and fourth column of Table. 
1. The entries in Table. 1 correspond to the choice of CKM parameters (p = 0.15, 77 = 0.33) 
and for 7/ in the final state we also include a second number in the parentheses corresponding 
to the second set of CKM parameters though the data seem to favor a positive value of p. 
The branching ratio for B — > rjK is suppressed by about a factor of 10 or more relative to 
B — ► vfK while B — > rjK* is enhanced relative to B — > rfK* by a small amount. This is 



in qualitative agreement with Ref [17] but we find a smaller enhancement B — > r/K* relative 
to B — > rj'K*. From the branching ratios in Table. 1, it is not possible to rule out the four 
quark operator explanation for the large branching ratio in B~ — > rfK~. 

For the quasi-inclusive decay we will use the same parameters as used in the exclusive 
decays with the second set of the CKM parameters as this set gives the larger rate between 
the two choices. We plot the decay distribution dT/dM rec in Fig. 1 showing the contributions 
from the effective Hamiltonian of four quark operators. From Fig. 1 we find the contribution 
from the four quark operators to the branching ratio B — > rj'X is around 1.3 x 10~ 4 for 
E, q i > 2.2 which is the signal region. This is far too small to account for the signal observed 
at high X s mass. This continues to be true even if when all parameters are allowed to vary 
over a reasonable range. 



5 Conclusion 



Summarizing, we have calculated the the effects of the effective Hamiltonian of four quark 
operators to the exclusive and quasi-inclusive decays of the B meson to 77'. Our analysis 
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indicate that the exclusive data can be explained by the four quark operators without sig- 
nificant contribution from the mechanism b — > sg*, g* — > f]'gg or from the intrinsic charm 
content of the rf . The contribution to the quasi- inclusive rate from the four quark operator 
is not enough to account for the observed signal. 
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Figure Caption 

Fig. 1 This figure shows the decay distribution as a function of the recoil mass M rec . 
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Branching Ratio |H| (BR) 
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/l) /8 
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(7.1±l;f d= 0.9) x 10~ 5 
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